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PACS 75 . 47 . Lx - Manganites 

PACS 71 . 27 . +a - Strongly correlated electron systems 
PACS 75 . 30 . Kz - Magnetic phase boundaries 

Abstract. - The field induced switching of conductance in the charge ordered half-doped man- 
ganites is controlled by the combination of metastability, an inhomogeneous high field state, and 
cation disorder. We study this non-equilibrium problem via real space Monte Carlo on a disor- 
dered strong coupling model appropriate to the manganites. We reproduce the variation of the 
switching fields with the mean ionic radius ta and cation disorder cta, and demonstrate how the 
experimental features arise from the proximity of several phases in the Landau free energy land- 
scape. Our prediction for the field melted state is consistent with a growing body of experimental 
evidence. 
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The manganites owe their fame to the 'colossal mag- 
netoresistance' effect, whereby an applied magnetic field 
strongly suppresses the resistance of the material. The 
most extreme version of this happens when the field can 
.ctually induce an insulator- metal transition (IMT). This 
occurs in the low bandwidth half doped manganites, e.g^ 
Pro.sCao.sMnOa, that are typically charge ordered (CO) 
insulators with CE magnetic order [1-3]. The CE-CO- 
I state is transformed to a ferromagnetic metal (FM-M) 
by an applied field through a first order phase transi- 
tion. The magnetic field induced conversion of the CE 
order to FM, and the crucial resistive switching, involve 
several subtleties which are only gradually being appre- 
ciated [4-14]. These include the surprising difference be- 
tween the melting fields in two lanthanide (Ln) families, 
namely the Ca [4] and the Sr [5,6] families; the smallness of 
the field melting energy scales [4-7] ; and the spatial char- 
acter of field induced melting and its relation to structural 
disorder [8-12]. We expand on these issues below. 

(i) The phase boundary between the CE and FM phase 
is first order, so metastable states exist on either side. 
This leads to hysteresis, and the switching fields are de- 
cided by non- equilibrium effects, not free energy balance, 
(ii) While trapping into metastable states is expected in a 
first order transition, an applied field leads to new equilib- 
rium states which are not simple continuation of the zero 
field state. In particular, it is increasingly apparent that 
the field induced conducting state is an inhomogeneous 



FM-M, a percolative metal with a finite volume fraction 
of charge ordered regions [8-12]. Finally, (iii) the field in- 
duced switching and inhomogeneous metallicity play out 
on a structurally disordered background. One expects the 
CO-I state to be more stable in narrow band systems. 
While this is true in the Ca based manganites [4], where 
the cation size mismatch cta is small, it is completely op- 
posite in the Sr based systems [5,6]. The more insulating 
system can be driven metallic easier in the presence of 
disorder! 

These three actors, metastability, field induced inho- 
mogeneity, and pinning disorder, enrich and complicate 
the field melting problem. There is no understanding of 
their interplay in the field induced IMT, and the non- 
equilibrium physics of the manganites. 

We study a coupled electron-spin-lattice model [15] ap- 
propriate to the manganites, in two dimensions (2D). For 
parameters where the FM-M and CE-CO-I phases are in 
proximity, we discover the following, (i) In the 'clean' 
limit, the forward (CE to FM) and reverse (FM to CE) 
switching fields, and h~ respectively, increase rapidly 
with decreasing bandwidth (BW), and below a threshold 
BW the field cannot melt the CO. This explains the results 
in the 'clean' Ca based manganites. (ii) The switching 
fields are dramatically affected by disorder: hf^Oas the 
bandwidth reduces, as in the Sr manganites, even though 
the CO is notionally stronger, (iii) The field melted state 
in the narrow band case is inhomogeneous even without 
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disorder - at finite field the system exists as a patchwork of 
metallic and CO regions, and the 'metallicity' is of percola- 
tive origin. Recent experiments have indeed verified coex- 
istence of FM-M and CO-I, e.g, in Lao.sCao.sMnOs [11] 
and Pro.sSro.sMnOa [9, 10], and the phase volumes are 
tunable depending on the field cooling protocol. 

We consider the field response of a model Hq that we 
had studied earlier [15] to capture phase competition. Our 
model is: 



af3 



Hn = 



{ij)a i (ij) 



(1) 



Ho involves two Cg levels per Mn, with nearest neighbour 
hopping [16], and Hunds coupling Jh between the Cg 
electrons and the t2g derived core spins. We assume a 2D 
square lattice. The Cg electrons have a coupling A to Jahn- 
Teller phonons, while the core spins have an antiferromag- 
netic (AF) superexchange J between them. The stiffness 
of the phonon modes is ji is the chemical potential, 
and the magnetic field, h = zh^ couples to the core spins. 
The 'phonons' are treated in the adiabatic limit and pro- 
vide a static distortion background for the electrons. The 
core spins are also treated as classical [17]. At half doping, 
n = 0.5, the clean model has a variety of phases of which 
we are interested in the CE-CO-I phase and its response to 
a magnetic field. Cation disorder leads to 'site disorder' 
due to the potential generated by the randomly located 
rare earth and alkaline earth ions, and 'bond disorder' in 
tij and J. We retain only the local random potential since 
it couples to the electron density and has a direct impact 
on charge order. It was recently shown [18] that while 
off-diagonal disorder is crucial for obtaining the spin glass 
phases at half doping, it is diagonal disorder that plays a 
dominant role in weakening and eventual disruption of the 
CO state. We model the potential seen at the Mn site 
as binary distribution P(e^) = ^{S{ei — A) + S{ei + A)) 
with variance A^. For Hq we use Jn/t oo, and set 
|S^| = 1. t = 1 is set as the reference energy scale and 
we also set K = 1. We restrict ourselves to J/t = 0.12. 
The chemical potential ja is adjusted to keep the electron 
density at n = 1/2 which is also x = 1 — n = 1/2. 

We use the travelling cluster (TCA) based Monte Carlo 
method [19] to solve the coupled electron-spin-phonon sys- 
tem above. Classical Monte Carlo techniques have been 
extensively used in literature [20] to study first order tran- 
sitions and the associated hysteresis. We will explore the 
dependence of the melting process on bandwidth (A/t), 
disorder (A/t), and temperature T/t. The experimen- 
tal variation of bandwidth (via ta) implies simultaneous 
change in A/t, J/t and A/t, while we have varied A/t, J/t 
and A/t independently. Our primary variation is in X/t 
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Fig. 1: Colour online: h — T phase diagrams, (a)-(b) Experi- 
mental phase diagrams in (a) moderate BW Lao.sCao.sMnOa 
and (b) narrow band Pro.sCao.sMnOs. Both materials have 
low cation disorder, a a 10~^A^. (c)-(d) Our results for two 
values of A/t and disorder A/t = 0. The parameters are chosen 
to mimic the experimental regimes above (see text) in terms 
of hysteretic response. Both FM-M and AF-CO-I states are 
locally stable in the shaded region. 



and, crudely, we compare this to experimental BW vari- 
ation. We have commented on this in more detail else- 
where [21]. 

All our results are for cooling at = followed by field 
sweep. The 'upper' and 'lower' switching fields, hf^ in the 
hysteretic response refer to the resistive transitions. 

The field induced melting of CO has been studied within 
mean field theory (MFT) earlier. For a model involving 
double exchange and Coulomb effects Mishra et al. [22] 
showed that suitable parameter choice could lead to a low 
he ~ 6T. Fratini et al. [23], considered a more elaborate 
model and worked out typical h — T phase diagrams within 
MFT, and Cepas et al. [24] used a variational approach at 
T = to estimate he- These studies, unfortunately, do 
not shed light on the spatial character of melting or the 
impact of disorder. 

Let us first clarify the field melting and hysteresis in 
manganites with weak cation disorder, cta is small, ~ 
10~^A^, in the lanthanide (Ln) family Lno.sCao.sMnOa [4] 
so as reference we look at Lao.sCao.sMnOs (LCMO), with 
a relatively large BW, and Pro.sCao.sMnOs (PCMO), with 
a smaller BW. Fig.l.(a)-(b) reproduces the h — T phase 
diagrams for these compounds from [4]. At low tempera- 
ture /i+ is - 20T for LCMO and - 30T for PCMO. The 
field h~ ^ at which CO is recovered on field reduction, is, 
however, < IT in LCMO, but - 15T in PCMO. h' 
at low temperature in LCMO indicates that the FM-M is 
metastable even at h = 0^ while at smaller BW (PCMO) 



p-2 



Field melting in CO manganites 




Nd Sm 




1.55 1.6 1.65 



1.55 1.6 1.65 
A/t 



Fig. 2: Colour online: Response to field cycling at A/t — 1.6 
and T/t = 0.02. (a)-(b) Clean system A/t = 0: (a) field re- 
sponse of the CO volume fraction (see text) and the FM struc- 
ture factor Sfm, (b) the resistivity p{h). (c)-(d) Corresponding 
results in the disordered case A/t = 0.12. Shaded region is the 
hysteretic window. Contrast the abrupt switching in the clean 
limit with the broad crossover in the disordered case. Note, 
both the CO volume fraction and the FM magnetic structure 
factors are normalized to unity. 



Fig. 3: Colour online: Tco, and hf at low T - comparing 
experiment and theory, (a)-(b) Data from experiments: (a) 
the Lno.sCao.sMnOs family, with typical a a 10"*^^^, (b) 
the Lno.sSro.sMnOs family with typical a a 10~^A^. Notice 
the rapid increase in hf with decreasing ta at 1-29 ^ in 

(a) , and the collapse of with decreasing ta at ta 1-34 in 

(b) . (c)-(d) Our results on the A/t dependence of Tco, Tce 
and hf. (c) Clean limit A/t = 0, (d) disordered systems, 
A/t = 0.12. The lines are a guide to the eye. 



the FM-M is no longer met ast able at = 0. The Tco 
at = for these two systems are ~ 220K and 250K 
respectively. 

Panels (c)-(d) show h — T phase diagrams computed 
by us at two couplings. For the lower value, A/t = 1.55, 
Tco/t ^ 0.05, while for A/t = 1.60, Tco/t ^ 0.06. Using 
a crude factor of 3/2 to convert from 2D to 3D and set- 
ting t ~ 0.25eV, these would correspond to about 210K 
and 250K. More significantly, A/t = 1.55, which is closer 
to the (FM-M CO-I) phase boundary, shows h~ = 0, akin 
to LCMO, while at A/t = 1.6, h' is finite and the FM- 
M state is no longer metastable at = 0. Although our 
scales, ~ O.lt, are larger than in experiments [25], the 
typical value ^{h'^ -\-h~) varies roughly similarly in the ex- 
periments and theory, changing from lOT to 20T between 
LCMO and PCMO and - 0.05t and 0.08t in Fig.l.(c)-(d). 
It is possible that including electron-electron interaction 
may reduce the melting fields somewhat [23] as mean field 
studies suggest. The h — T phase diagram is derived from 
the evolution of the q = {0, 0} feature in the magnetic 
structure factor 6'(q) and the 'volume fraction' Vco [26] of 
the CO phase. 

Fig. 2 shows these features at A/t = 1.6 and T /t = 0.02. 
Fig. 2. (a) shows the 'switching' of Vco and Sfm = 5'(0, 0) 
on field cycling, for a clean system. The complementary 
switching in FM and CO volume fraction is sharp and si- 
multaneous, along with the transition in resistivity p{h) 



shown in Fig.2.(b). Notice the residual Vco in the field 
melted metallic state in Fig. 2. (a). For a weakly disordered 
CO state, at A/t = 0.12, Fig.2.(c) shows that the CO 
melts on increasing /i, but fails to recover on downward 
sweep. The magnetism in Fig.2.(c) and p{h) in Fig.2.(d) 
still show hysteresis but the transitions are quite broad. 
This suggests that the field induced 'melting' in the dis- 
ordered case occurs gradually^ unlike the abrupt switching 
in the clean case. 

Using indicators as shown in Fig.l and Fig. 2 we made a 
comprehensive comparison between experiment and the- 
ory for the dependence of transition temperatures and low 
T melting fields on BW and disorder. Fig. 3. (a) is data 
from the 'low disorder' Lno.sCao.sMnOs [4] family, and 
Fig.3.(b), is for the moderate disorder Lno.sSro.sMnOs [5] 
family. We compare these with our results on the A/t de- 
pendence of Tco and in Fig.3.(c) and Fig.3.(d), with 
A/t = and A/t = 0.12 respectively. 

(A) In the clean limit. Fig. 3. (a), Tco increases with de- 
creasing BW, Tce decreases, and the fields increase 
rapidly beyond Smo.sCao.sMnOs (SCMO). A comparison 
of Fig. 3. (a) with Fig.3.(c) shows that our results capture 
all the experimental trends in the clean limit. The increase 
in Tco is due to the decrease in kinetic energy with de- 
creasing rA (increasing A/t in our case), while Tce is sup- 
pressed because CE magnetic order is also kinetic energy 
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Fig. 4: Colour online: Spatial snapshots of local charge order 
(top row) and nearest neighbour magnetic correlation Si.Sj 
(bottom row) on field cycling, X/t = 1.55, A/t = 0, T/t = 0.02. 
/i = (left) initial state, h > ht (center), and /i = (right) 
after the field sweep. The lattice size is 40 x 40. In the bottom 
panel the red bonds are FM the green bonds are AF. In the 
top panel grey regions imply metallic regions and the red and 
green regions are the 0.5 — ^ and 0.5 + ^ onsite charge densities. 



driven and gives way to a 'G type' {7r,7r} phase at large 
A/t. The increasing separation of Tco and Tce implies 
that the charge order becomes progressively independent 
of CE magnetic order with decreasing BW, and so less 
responsive to an applied field. This is precisely what one 
observes for in Fig. 3. (a) & (c), where below a critical 
BW the CO no longer melts! In fact for \/t > 1.60 we 
observe that the high field state is FM-CO and not FM- 
M. The magnetic phase change no longer drives a resistive 
transition. 

(B) For the disordered family, Fig.3.(b), Tco and 
decrease rapidly as we head towards small BW 
Smo.sSro.sMnOs (SSMO). Our results in Fig.3.(d) probe 
BW variation at disorder A/t = 0.12, and capture the 
correct dependence. The fall in Tco at large A/t is due 
to the competition between the 'stiffness' {Keff{X), say) 
of the CO state, which weakens at large A/t [27], and the 
pinning effect of the random potential e^. This is similar 
to the random field Ising model [28] where the exchange 
(Jeff) serves the role of our 'stiffness' and the random field 
{hi) is akin to e^. As is well known, long range order is 
lost below a threshold for Jg/ / / , equivalent to our 
K(.ff/A. The weakening and ultimate destruction of the 
zero field CO state with decreasing (or t/A) automati- 
cally suppresses the melting fields hf . The non monotonic 
character of Tco and hc^ with a reduction also at the large 
BW end, is due to the presence of the competing FM-M 
phase at A/t < 1.5. 

While the trends are well captured, let us compare the 
experimental Tco scales with theory. At the 'large' BW 
end, rA ^ IMA, 'clean' LCMO has Tco ^ 220K, while 



Fig. 5: Colour online: Landau landscape. We plot the free 
energy F((/), h), where is the 'order parameter' (on the x axis) 
and can have FM-M, CE-CO-I, AF-M, or AF-CO character, h 
increases from hi to /is (on the y axis). We have not shown the 
metastable FM-CO state, for /i ^ /i2, to avoid cluttering. 



'disordered' Pro.sSro.sMnOa (PSMO) has a suppressed 
Tco ^ 150K. If we identify A/t ~ 1.55 crudely as the 
equivalent, the clean Tco is ~ 0.055t ~ 225K, and the 
'disordered' Tco is ~ 0.035t ~ 145K. The peak Tco is 
- 160K in Fig.3.(b) and it is - 0.05t - 200K in theory. 
Given that we are solving a 2D model and have not fine 
tuned parameters the correspondence is reasonable. 

The thermodynamic indicators, however, do not reveal 
the spatial character of the 'melted' state, while Fig. 2. (a) 
& (c) suggest something unusual: there is finite CO phase 
fraction in the 'FM-M' ! Our map of the density field in 
Fig. 4 reveals how the system is spatially organised. The 
lower panels show that the CE phase transforms to a FM 
at large field (and only partially recovers CE order on re- 
verse sweep) but the large h density field is strongly inho- 
mogeneous. The middle panel, top row, shows charge mod- 
ulated 'minority' regions in a homogeneous background, 
surviving way past the CE FM transition. For the pa- 
rameter values in Fig. 4, and overall for A/t < 1.6, the 
'high field' state is a percolative metal. Beyond A/t = 1.6 
it is an insulating FM-CO. 

The peculiar high field state requires us to examine the 
various phases that 'compete' in the manganites around 
X - 0.5 at finite h. At = the x = 0.5 CE-CO-I is 
separated from a FM-M phase at lower x, and an AF-M 
phase at larger x, by windows of PS [15]. However, beyond 
a small threshold field a homogeneous x = 0.5 state is no 
longer possible [21] and the x = 0.5 system breaks up into 
either (a) x < 0.5 FM-M and x > 0.5 AF-M patches, 
for A/t < 1.6, or (b) x < 0.5 FM-CO and x > 0.5 AF- 
M patches for A/t > 1.6. These intermediate field states 
evolve, respectively, into a homogeneous FM-M or FM- 
CO state at large h. While the finite h equilibrium state 
is complicated by PS, field sweep experiments (and our 
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simulation) also probe the existence of metastable states 
in which the system can get trapped. Fig. 5 summarises 
the situation at \/t ~ 1.55 in terms of a Landau free 
energy schematic. At /i = the CE-CO-I is the absolute 
minimum while the FM-M is metastable (ref Fig.l.(c)). 
This continues to finite /i, we show a typical field h = 
hi. At h = /12(A) the x = 0.5 state phase separates into 
AF-M and FM-M, and there is also a metastable FM- 
CO. Finally, at some h = /13(A) the equilibrium state is 
a homogeneous FM-M but the continuing presence of the 
low energy metastable FM-CO state would affect all low 
T field sweep experiments. 

To conclude, we have clarified how non-equilibrium ef- 
fects and disorder control the field driven conductance 
switching and insulator-metal transition in the half doped 
manganites. We discover that the field induced conduct- 
ing state is spatially inhomogeneous, even in the absence of 
quenched disorder. Our results on the h—T phase diagram 
and the spatial character of the melted state address ex- 
isting results and growing body of spatially resolved data. 

We acknowledge use of the Beowulf cluster at HRI. 
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